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S’aiOPSIS 


The present work is devoted to the investigation of the 
linear stability of inviscid parallel shear flows* 

Howard (1961)/ in one of the inost beaut if 'ol theorems in 
hydrodynamic stability of homogeneous shear flows/ has shown 
that the eigenvalues c of unstable modes lie inside a semi- 
circle/ in the upper half of the conplex plana, whose diameter 
is the range of the basic velocity* This result does not 
incorporate the effect of curvature of the basic velocity 
profile on the stability of the homogeneous ^ear flows* But, 
it is curvature of the basic velocity profile that plays 
dominant role on the stability or instability of homogeneous 

I 

shear flows (Rayleigh, 1890; Pjortoft, 1950; Tollmien, 1935, 
Lin, 1955)* Motivated by these arguments/ the problCTi is 
re-investigated and an instability region depending on the 
curvature is obtained* This region comes out to be dqpth 
dependent as well* In conjunction with Howard's result, it 
reduces the instability region given by Howard's semicircle* 

Instead of a homogeneous fluid, if one considers a 
horizontally stratified fluid under the action of gravity, then 
the question. What is the role of curvature on the stability 
of the flow? does not have a satisfactory answer till today* 

In fact/ the validity or otherwise of the Rayleigh's inflexion 
point theoirem in the presence of stratification has not been 



established so far* Motivated by these argximents, the role of 
curvature of the basic flow profile on the stability of 
stratified shear flows is studied. Simple general i25at ions 
of Hoiland^s (2 951) criterion and Sattinger's estimate are 
given- In addition to generalizing the instability region 
discussed above/ it is shown that the complex wave velocity 
of an arbitrary arbitrary unstable mode lie inside a samiellips* 
region whose major axis coincides with that of Kochar-Jain'' s 
(1979) semiellipse/ while its minor axis d^ends not only on ! 
the stratification but also on the curvature of the basic I 

velocity profile# 

1 

The work of Howard (1961) show that the conplex wave i 

velocities of unstable and their conjugate damped nK>des lie 
inside a bounded region in the Cj^c^-plane# The wave velocities 
of singular neutral modeS/ by their very definition, lie between 

the miniimm and maximum values of the basic velocity# In this 

' - ^ 

work/ it is shown that the wave velocities of non— singular 
neutral modes of the Tayior-Goldstein problem are also botmded* 

Stability of parallel shear flow of a compressible fluid 
is of interest in meteorology# Unfortunately however, there is 
not much systematic work on the stability characteristics of 
inviscid compressible parallel shear flows# Following Blxomen 
(1970) and Chirwonas (1970), this problem is studied with and 
without the presence of gravity# Following Blrnien gravity is 
ignored and the basic thermodynamic state is assuiited to be 



constants It is shown that shear free basic flow (U S o) 
sv 5 >ports supersonic waves- Further^ it is shown that the 
corrplex wave velocity of an unstable subsonic disturbance 
lie inside a semiellipse type region depending on the Mach 
number M* Eckart^s (196 3 ) semicircle comes out to be a 
special case of this when M = 0» Incidently/ this region 
comes out to be depth as well as wave number d^endant» For ; 

I 

unboiinded flows^ a sufficient condition for stability to 
supersonic disturbances and an estimate for the growth rate | 

I 

I 

of an xinstable supersonic disturbances are also aiven* i 

i 

I 

For stratified conpressible shear flows/ it is shown that 
the instability region for subsonic disturbances is an 
semiellipse type region, which depends on the Richardson number,; 

wave number and depth of the fluid layer* For flows with 

# 1 '' 
^min ^ this region reduces to the line c^ = O when '^0 ■* 4 *“ 

in accord with Miles' theoron (Chimonas - 1970 )• Under an 

approximation due to Shivaraoggl (1977), the role of curvatxire 

on the stability is also studied* 



CHAPTER 1 


introduction - 

St ability Probi Am ^ 

ConsMer the notion of an inviscid oo^^esslhle fluid 
tonfined between two horizontal infinite rigid planes. The 
governing equations are the Euler equation 


^3t 

the equation of 


9P 

^ + (U*v)p + pv.U _ o. 


^ (H + (u-v)u) = - vp , 


( 1 ) 


continuity 


( 2 ) 


and the equation of state 
P = pCp^s) 

(3) 

^or sinplicit., we shall consider only isentropio flows fi,r 
Which 


3 S . a 

^ + (U*v)s S o^ 


(4) 


where 3 is the velocity, p the density, p the 
s the entrc^y of the fluid and g is the 


pressxire - 


gravity. 


acceleration due to 


Rations (3), (4) end ( 2 ) can be ocnfcl^^j 

ap 9 

al + + aj s= o 


( 5 ) 
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2 3 P 

where a,^ = is the square of sound speed* The bovtndary 

conditions are that the vertical conponent of the velocity- 
vanishes on the rigid walls at y = yi^y 2 * 

If we fix a coordinate system with the x— axis along the 
horizontal direction and y-axis along the vertical direction, 
then we see that the flow variables given by 

U = U(y) i, p = p (y) and p ss p(y), (6) 

where i is the unit vector along the x»-axis, satisfy the 
governing equations (l)# (2) and (5) and also the boundary 
conditions provided 

^ - pg (7) 

Here, we allow U(y), p(y) and p(y) to be any twice continuously 
differentiable functions* These fields define the basic flow* 
''Yet not every solution of -the equations of motion, ev^ if 
it is exact, can actually occur in Nature* The flows that 
occur in Nature must not only obey -the equations of flxiid 
dynamics^ but also be STABLE* Landau and Lifschltz (1959)* 

If the basic flow is disturbed slightly, the disturbance 
may either die away, persist as a disturbance of similar 
magnitude or grow so muc^ that the basic flow becomes a different^ 

i 

laminar or a turbul^t Oow* Broadly peaking, we call such 1 
disturbances { asyc^toticaliy) stable, neutrally stable or 
unstable respectively* A3JL po^ible slight disturbances are ' 
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likely to be excited in some degree by small irregxxlarities 
or vibrations of the basic flow in practice, so the basic 
flow will persist only if it is stable to all slight 
disturbances# Mathematically, we can define stability in the 
sense of Lyap\inov, but in most applications conanon sense makes 
a formal definition unnecessary# 

Now, we shall study the stability of the basic flow 
given by (6) and (?) to two-dimensional, in finitesimal, wavy 
disturbances# Such disturbances are called normal modes# 

Let the dist-urbed motion be given by 
IT = (U(y) + u^(x,y,t), v'Cx,y,t)) 

( 8 ) 

P as p(y) + p^(x,y,t) , p = p(y)+p^(x,y,t) 

For infinitesimal dist\irbances the equations can be 
linearized and the linearized perturbation equations are 

-I# ‘9) 

+ =^||r-p'g (10) 

For wavy disturbances we can take the normal modes 



4 


u' 

II 

exp 

[ik (x 

- ct)] 




v' 

= v(y) 

exp 

[ik (x 

- ct)] 







[ik (x 

ct) ] 



(13) 

p' 

= P(y) 

exp 




p' 

= p(y) 

exp 

[ik (x 

*• ct)] 




equations (9 

) - 

(12) lead to the 

equations 


p [ik (U - 

c)u 

+ U'v] 

= - ik p 



(14) 

p [ik (U - 

c)v 

] =-P^ 

^ - pg 



(15) 

ik 

(U ^ C 

)p t 

P' V + 

pCik u + 


* 0 

(16) 

ik 

(U - c 

) p - 

ry 

- pg V + pa^ [ik 

/V 

U 4* 

O'] = 0 

(17) 
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is the Brunt-Valssala frequency# 
Equation (18) can be rewritten as 


_r(U-c)H^ - rU^H , 
[ •; • ] 


- rk"^ (Ui-c) H + 


rN^H 

Tu-cJ 


= 0 


( 20 ) 


where H(y) = v(y) exp (- j ) 


y 2 q 

and r(y) = (y) exp ( / dy') 


(21) 


The associated bovindary conditions are 


H(y^) = 0 = HiY 2 ^ 


( 22 ) 


The phase velocity o = (c^ + Ic^) „ay oon*,lex and the 
wave number k is taken to be positive, it is assumed that 

N^y) > o vSiloh itKans that the stratification is statically 
stable# 


The stability equation ( 20 ) and equation (22) define 
an eigenvalue problem for the phase velocity c. If all the 
eigenvalues are real, th^ the basic flow is stable to two- 
dimensional, infinitesimal, wavy disturbances# The flow is 

unstable if there is a ccHi^lex eigenvalue c s= (c + ic ) with 

,X' ' ^ 'i ■ ■■ 

c jL O ♦ 

Let H be an unstable solution. Then the transformation 
H e (u-c)f transforms (20) to 



6 


^r(U-c)^ ^ 2 v2 2 

[— - rk^ (u-c) P + rW^P = 0 


(U-e) 

5- 


(23) 




and the associated boundary conditions are 


P(y. ) = O = PCYo) 


(24) 


Since c^ > O for an xinstable mode* a proper meaning can be 

attached to tine transformation 

1 

G = (U-c)^P (25) 

by taking a definite branch of it* The transformed equation 
in terms of G is 


r(U-c)G* tr 

[——53 - I [ 

. (U-c)? ^ 

1 ^ 


rU^ 


^ (U-c) 

1 *• 5 - 


j] G - 


rU'^G 


4(U-c)(l - 


(u«-c)- 


- rk2(u-c)G t ^ = 


0 


(.26 ) 


and the associated boundary conditions are 


G(yj) = O = G(y2) 


(27) 


2* Homogeneotis Shear Plows 

For incon^ressible constant density fluids the stability 
equation is tine Rayleigh equation 


H» - k:^ - 




H 


(28) 
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and the boxindary conditions ajre the same as (22)# In 1880^ 
Rayleigh proved his famous inflexion point theorem namely^ a 
necessary condition for instability is that vanishes 
atleast once in the flow domain* This result was extended 
by Fjortoft (1950) and Holland (1951 )* None of these conditions 
are sufficient conditions* as is evident from the Gountei:exait|>le 
of Tollmien (1935)» However* Tollmlen (1935) and Lin (1955) 
have shown that these conditions are sufficient for the 
instability of symmetric profiles* Howard (1961) has proved 
the semicircle theorem by \^ich the conplex wave velocity for 
any unstable mode must lie inside the semicircle in the Tapper 
half plane with the range of basic velocity as diameter* 

Holland (195l) has also found an estimate for the growth rate 
of an tmstable- Sattinger (196?) has fotind an estimate for 
the growth rate* which depends on the curvatxire of the basic 
flow profile and d^th of the fluid layer# In the stability 
of homogeneous shear flows* the curvature of the basic -v^ocity 
profile and the depth of the fluid layer play important roles* 
For a thorough discussion of this problem* one may be referred 
to Lin (1955)* Drazin and Howard (1966) and Drazin and 
Reid (1981)# 

3# Stratified Tnconuressible Shear Flows 

The stability of i^ear flows of a continuously str^ified 
fluid is a fascinating phenomoion of great iitportance to 
meteorology# Helmholtz (1968) was the first to study the 
instability of an interface dividing two layers of invi^ld 
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and Incortpressible fluids of different densities and in 
relative horizontal motion# However, the detailed and 
exhaustive study of this problen was done by Kelvin (1871) 
in context with the study of the generation of water waves* 

The stability equation for this problem is 

(PH' )' - e A - H + = 0 (29) 

which is obtained from ( 20 ) by taking = 0 in it* The above 

equation is known as the Taylor*-Gold stein equation in honour of 
their derivations and exploitations in i931« 

The first paper with general velocity and density 

distributions was published by Synge (1933). He generalized 

Rayleigh's theorem and obtained bounds for the growth rate 

which were later obtained Independently by Yih (1957) and 

Drazin (1958). Miles (1961) proved a conjecture of Taylor 

namely the flow is stable if the Richardson number J > 4 

— 4 

everywhea^e in the flow domain for a class of flows# Howard 

I 

(1961) gave a sinple proof of Miles' theorem aixi removed the 
restrictions on the velocity and density profiles# He alaa 
proved that for any xinstable mode,^ the complex wave velocity 
must be inside tte semicircle in the r^sperhalf plane, which 
has the range of basic velocity as diameter# Kochar and Jain 
( 1979 ) improved upon this result by showing that the complex 
wave velocity for any unstable mode lies inside a sCTiielllpse 
rather than a semicircle whose major diameter coincides with 
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Howard^ s semicircle^ while minor diameter depends on stratifi- 
cation* This result was further generalized by Jain and 
Kochar (1983) and Makov and Stepanyants (1984) incorporating 
the wave mamber and the depth of the fluid layer# Drazin 
(1958) found an exairple which is stable when J > 4 and xanstable 
when *1 4 * 

Howard (1963) generalized the method of perturbing the 
stability boxandary and finding the adjacent xinstable modes/ 
to the case of stratified flow* Hupp e3rt{ 1973) considered 
number of exanples and established the validity of Howard's 
formula* Engevik (1973/ 75/ 78) gave a sirtple derivation of 
Howard's formula and derived the differential equation to find 
the stability boiindary* Engevik et al# (1985) have discussed 
the perturbation of the stability even when Howard's formula 
is not aj^licable# 

After the early work of Synge (1933)/ many workers have 
discussed the role of curvature of the basic velocity profile 
on the stability of stratified shear flows# Thorpe { 1969 ) 
found a counterexanple for Ejortoft's theoroDO* Yih (1974) 
gave his sufficient condition for stability/ viiich is the 
strongest known generalization of Raylei^'s theorem# Saner jee 
et als* (1972/ 74/ 78) obtained instability conditions depending 
on the curvature of the basic velocity paofile for certain class 
of profiles# 
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Banks et al» (1976) examined the overall pattern of 
normal modes. For a given flow and wave number the modes 
are divided into five classes some of which may be ^ipty* 

4- Compressible Shear Flows 

Stability of parallel shear flow of a conpressible 
fluid is of interest in meteorology* Landau (1944), Hatanaka 
(1947) and Miles (1968) studied the instability of a vortex 
sheet in a gas with respect to infinitesimal disturbances 
neglecting all diffusive processes* But Bltimen et al* (1975) 
have casted dotibt on the physical value of models ^ich 
incorporate a vortex sheet* Unfcrtunately, however, there is 
not much systematic work on the stability characteristics of 
inviscid coxrpressible parallel shear flows (Betchov and 
Criminal e, 1967, Dandapat and Gupta, 1977) • 

Some general stability characteristics of ga^ynamic 
shear flows, without including the external force field of 
gravity, have been ortained by Lees and Lin (1946)* According 
to them, the distTirbances can be classified as subsonic, sonic 
or supersonic d^o^ding on vdiether the relative phase ^eed 
lu-cl is less than, equal to or greater ths® the sound sp®^i# 
i.e. lu-cl ^ a^* The siibsonic disturbances are the coi3nte2>- 

parts of the so-call^ inertial modes# which are elutions 
of the corresponding equations for inconpressible fluids* The 

i • 

physical significance of the sonic disturbances is ^parently 
not dear* Finally, si: 5 )ersonic disturbances correspond to 
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coiipression or sound waves* Eckart (1963) generalized Howard's 

semicircle theorem for gasdynamic shear flows# Chimonas (1970) 

extended Miles' theorem and Howard's estimate for the growth 

rate of unstable modes to gasdynamic shear flows* The extension 

of the Rayleigh stability criterion and Howard's semicircle 

theorem to coirpressible shear flows^ obtaincad by Lees and 

Lin (1946) and Eckart (196 3) respectively, have each been 

rederived by Blximen (1970) by a different approach* He also 

found a subsonic neutral solution of the stability equation 

when the basic flow is represented by the hyperbolic'-tangent 

velocity profile and discussed xmstable modes ntamerically. 

Shivamoggi (1977) studied this problsn under the spproxiraation 
i 

— ~ « 1 and Cj « 1* The works of Chimonas and Bivins have 
a^ i 

been extended to Spiral flows by Howard (1973) and Dandapat 
and Gupta (1975) • 

5* General Aspects 

Linear stability analysis can indicate only instantaneous 
tendency of a laminar flow for small pertvirbations* The 
noment the flow is unstable and ^ows a t^adoacy to grow the 
linear analysis will no moire be valid* Sizable perturbations 
and Instabilities subsequ^it to primary instability become 
inportant* But linearized analysis must be considered as the 
first stqp in any stability theory and it is a natural 
starting point for the description and definition of 
nonlinear problems* An Idaatioai relationship exists also 



12 


between stability of inviscid flows and viscous flows# in 
which the inviscid analysis is considered as the first step 
in any stability theory and natural starting point for the 
description and definition of viscous problems* For a 
discussion of these points one may be referred to Lin (1955)# 
Chandrasekhar (1961 ). Eckhaus (1965), Joseph (1976) and 
Drazin and Reid (19 81). Itirther, the instability mechanism 
is related to the over-reflexion of waves and this has beai 
discussed in Acheson (1976) and Pellacani (1983)* 

6 • Plan and Contributions of the Thesis . 

The present work is devoted to the investigation of 
the teirporai stability of inviscid parallel shear flows 
using the normal mode technique* 

In the seccxid chapter of this thesis we consider the 
reduction of Howard^ s semicircle bound on the range of the 
cortplex wave velocity of an arbitrary tmstable mode in the 
stability problem of homogeneous shear flows* As pointed 
out earlier, Hbwai^d (1961) has shovm that for instability 
of a homogmeous ^ear flow the cxjnplex wave velocity c of 
any arbitra 3 -y xmstable mode must lie inside a semicircle in 
the upper half plane, with the range of the basic velocity 
as diameter* This elegant result does not contain any 
parameter that characterizes the curvature of the basic 
profile or the depth of the fluid layer* The famous Inflexion" 
point theorm of Rayleigh Cl880) and its subsequent extensions 
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by Hoi and (1951) and Fjortoft (1950) show the in|>ortant TOle 
played by the curvature of the basic velocity profile on the 
stability of homogeneous shear flows* The estimates for 
growth rate of an unstable mode of Sattinger (1967) and Craik 
(1972) and also the exanple of Tollmien (1935) show the role 
of the depth of the fluid layer on the stability of homogeneous 
shear flows* In the present work Howard's semicircle botmd 
is fTorther reduced* The reduction depends on the curvature 
of the basic velocity profile and also on the d<^>th of the 
fluid layer* 

If, Instead of a homogeneous shear flow^ one considers 
a hetrogeneous shear flow^ then gravity plays a dominant role 
on stability characteristics. For this p3X>blenv the role of 
curvature of the basic velocity profile on the stability 
characteristics is not clear* In fact the validity of the 
res\ilts available for homogoieous shear flows, in the presence 
of stratification have not been established so far. This 
iTK>tivates us to study, in the ■yiird chapter, the role of 
curvature of the basic velocity profile on the stability of 
stratified shear flows. 

First, we give a simple generalization of Holland^ s 
result to the stratified case* Then, we generalize Sattinger' s 
estimate for the growth rate of an xinstable mode to the strati- 
fied case. This estimate Involves the curvature and we note 
that its ±ii^rovCT«nt by Craik (1972) cannot be carried over to 
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the stratified case. Recently Kochar and Jain (1979) have 
proved that the \instable nodes of a statically stratified 
shear flow lies inside a certain semiellipse in the upper 
half-plane^ whose major axis coincides with the diameter of 
Howard's semicircle, while its minor axis d^ends on the 
stratification parameter viz. the Richardson number. In oxir 
work, we show that the unstable modes lie also inside a 
semiellipse region whose major axis coincides with the 
diameter of Howard's semicircle, while its minor axis d^ends 
not only on the stratification tut also on the cuirvature 
of the basic velocity profile* Lastly, we extend the restilts 
of the previous chapter to the stratified case. It is 
interesting to note that these resiilts are also inprovements 
over the results of Saner jee et al. (1974, 78). 

In the fourth chapter, we discuss the boxandedness of tte 
phase velocities of the normal modes of the TaylorKSoldstein 
problem* The works of Howard (1961) and Kochar and Jain (1979) 
show that the conplex wave velocity of unstable and their 
conjugate damped modes lie in a bounded region in the 
plane. The wave velocity of singular neutral modes, by their 
very definition, lie betwe^ the minimum and maximum values 
of the basic velocity. But,v^at about the boundedness of the 
phase velocities of non-singular neutral modes? We show that 
the wave velocity of an arbitrary non-singular neutral mode 
of the Taylor *<301 dstein problm is also bounded. We give 
explicit bounds. 
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In -the fifth chapter, we study the stability of conpre- 
ssible shear flows without the presence of gravity# We show 
that the basic flow U S O supports supersonic waves but not 
subsonic waves* The semicircle theorem has he&i proved for 
compressible shear flows by Eckart (1962) and later by Bltimen 
(1970). But this resiilt does not incorporate the effect of 
compressibility and remains the same as that for inconpresslKLe 
fluids# We show that the complex wave velocities of an 
unstable subsonic disturbances lie inside a semiellipse type 
region depending on the Mach nTimber# The semicircle comes out 
to be a special case of this when M = O# Incidently, this 
region comes out to be depth as well as wave number d^endant* 
When the flow domain is unbounded, we have rederived Hoiland^s 
estimate for the growth rate of any unstable mode by a 
different method and have also found a sufficient condition 
for stability to supersonic disturbances# We show by an 
example that the phase velocities of non— singxilar neutral 
miodes of compressible flows are not always bounded# 

In the sixth and final chapter# we study the stability 
of parallel coc^ressible shear flows in the presence of gravity# 
We modify Chimonas' (1970) proof of Miles' theorem# This 
enables us to show that the instability region for subsonic 
disturbances is not Eckart' s ^micircle, but rather a semi- 
ellipse type region vdiich depends on the Richardson nuinber 
and also cmi the depth of the fluid layer and wave nuinber# 
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Furthermore# when ^ this region reduces to the line 

1 

0 ^ = 0 when "Iq 4 — es expected from Miles* theorem^ Following 
Shivamoggi (1977)# we take and the imaginary part of conplex 
wave velocity Cj|^ to be very small compared to unity so that 
their product can be neglected in coitparison to unity* Under 
this approximation# we study the role of curvatTore of the 
basic velocity profile on the stability of the flows* For 
subsonic disturbances, we have generalized many standard res^alts 
known for incoitpressible flows like instability criterion and 
estimate for the growth rate# due to Synge* 



CHAPTER 2 


ON REDUCING HOWARD^ S SEMICIRCLE FOR 
HOMOGENEOUS SHEAR FLOWS 


1# Introduction 

In the stability problem of homogeneous shear flows 
Howard (1961) established that the coitplex wave velocity 
c = (Cj. + iCj^) of an arbitrary \instable mode (c^ > o) 
must be inside or on the semicircle 

Cc^ - (1) 

where a =s and b = ^ the basic velocity profile# 

This elegant result does not aj^ear to have been iitproved in 
the literature. The famous inflexion point theorem of Rayleigh 
and its subsequent extensions by Pjortoft (1950) and Hoiland 
(1951) show the inport ant role played by the curvature of the 
basic velocity profile on the stability or instability of 
homogeneous shear flows. However, Howard^ s senicirde theorem 
does not contain any parait^er that characterizes the cxirvature 
of the basic velocity profile. In this chapter, it is shown 
that Howard's i^mlcircle bortnd on the conplex wave velocity 
of an arbitarary unstable mode can further be reduced# The 
reduction depends on the curvature of the basic velocity 
profile and also on the depth of the fluid layer. The role 
of the depth of the fluid layer on the stability of homog^aeous 
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shear flows have been brought out in the works of Sattinger 
(1967), Craik (1972) and Tollmien (1935). 

2 • Mathematical Analysis 

The basic equation and the boundary conditions for 
the problem of homogeneous shear flows are given by (cf» 
Chapter l) 

[(U-c)G']'^[^" + k^(U-c) + G = 0, (2) 


G(y^) = 0 = G(y 2 ) (3) 

where the primes denote diffeiientiation w«r»t* y, k is the 
wave number and c = c^ + ic^^ is the complex wave velocity* 


Theorem 1 * 


If (c,g)^ c = c + ic 


( 2 ) - (3) and f =- + 

then 


2b 71 


i’^ 

2 


(yi-Yj)^ 


c < X (c •• — ^) 
^ '‘^r m+1^ 


where X = [ Imax ^ 


> 0 is a solution of equations 
■ > O, for eveiry y e [y^^y^]. 


(4) 



a ^ O* 


Proof - 

Multiplying (2) by G* (the conplsx conjugate of g) and 
integrating over the range of jg we have \:^on integration by 
parts once and making use of (3)# 
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; (U-c) [ 1 G' j 1G|^] + |; U" IGI^ dy + / = O 

(5) 

The limits of integration and the Infinitesimal length dy is 
dropped throughout for convenience* Equating the real and 
imaginary parts of (5) to zero and cancelling c^^ (> O), we get 


/ (U-c^) [iG'i^* + k^lGf^] + <1 / U*' IGI^' 



U'^(U~c^) IGI^ 
[(U-c^)2 4- cj] 


(6 ) 


j [lon^ -f. - i : - g L-J iS i L — , o. (7) 

^ + oj] 

Multiplying (7) by mc^ and adding the resulting equation to (6), 
we get 


/ [ut(m-l)Cj^] [IG'I^ + k^lGl^] + I / U» IGI^ 


^ U'^ [U - (mfl)Cj.]lGl 

^ [{U-Gj.)^ + c|] 


= 0 


( 8 ) 


Since a > O, and a < < b (Rayleigh, 1330 ), it follows that 


[u+(m*-l = b > O and = -a < O (9) 

Equation (8) upon using (9) and the Rayieigh-Ritz inequality 
(Shultz, 1973) 


; IG'I^ > 


71 


Cy^-Yi) 


^ / IGI 


do) 



20 


gives 


[2cJ f - u'^ I (ml)c^ - y }] IGl^ 
4 [(U'x:^)\c^] 


< 0 


( 11 ) 


Under the conditions of the theorem, (ll) clearly Irplies 


that 




where 


. rU'^'- (rcrt-l)-, 

X = L 2f -I 


max 


and this establishes the th^rem* 


Theorem- 2# 


Under the conditions of Theorem 1, if 




mf-1 


( 12 ) 


then the parabola 




(13) 


intersects Howard's s«nicircle (l)» 

Proof * 

It is easily seen that the parabola (Is) touches 
Howard's semicircle (l), if 




ItH-l 


( 14 ) 
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The value of X given by (14) with the positive sign is 
rejected as it leads to < a which violates a < c^ < b# 
Hence^ if 

X < a+l.) - [(^ a+b)^- - (b-a)‘]V2 

then the parabola (I 3 ) intersects Howard's semicircle (l)# 


This proves the theorem* 

Theorem 3* 

If (c,G), c = c^ + ic^jr c^ > O is a solution of 
equations (2) - ( 3 ) and every y e then 


c^ < 


(c^ 



(15) 


where X* = [ max ^nd m = | , a > 0, 

Proof * Multiplying (?) by - adding the resulting equation 

to (6), and pixjceeding as in Theorem 1^ we get the result* 

Theorem 4 ^ 

Under the conditions of Theorem 3^ if 


X* < t*a) -[(^ t*a)^ - (b-a)^]l^^ (16) 

then the parabola 

4 = 


intersects Itoward's semicircle (i) • 
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Proof * Follows by proceeding as in Theorem 2* 

3* Concluding Remarks 

Theorems 1 -- 4 in conjunction with Howard's semicircle 
clearly show that if 

a > O and either f > 0 or g < (IS) 

for every y e [y^/y 2 l, then tinder the conditions (12) or (16), 
Howard's semicircle bound for honxjgeneous shear flows can 
further be reKluced# Further since the conditions (18) allow 
the basic velocity profile to change the sign of its curvattire 
somewhere in the flow domain the result assiames more significa- 
nce because under these conditions the basic velocity profile 
if symmetric is unstable (Tollmien, 1935)» 



CHAPTER - 3 


ON THE ROLE OF CURVATURE IN THE STABILITY 
OF HETROGENEOUS SHEAR FLOWS 

t » Introduction , 

For a steady, plane, parallel flow of an inviscid, 
incorrpressible, hortK>geneous fluid confined between two parallel 
rigid planes, a necessary condition for instability is that 
the basic velocity pirofila has atleast one inflexion point in 
the flow domain* This is the fcwnous Rayleigh's theorem* 
Hoiland (1951 ) extended this result further by proving that 
a necessary condition for instability is that U*(U-c^) < O 
atleast once in the flow domain, where U(y) is the basic 
velocity and c = c^ + ic^, > 0 is the phase velocity of 
an unstable normal mode* Fjortoft (1950) gave a stronger 
result by showing that a necessary condition for instability 
is U*(U — Ng) < 0 atleast once in the flow domain where 
Ug = U(yg) such that U*(yg) =* O* For a discussion of these 
results and their physical interpretation, one may be referred 
to Drazin and Reid (1981)* These results Illustrate the 
iirportant role played by the curvature of the basic velocity 
in the stability of homogeneous shear flows* 

However, instead of a homogeneous fluid, if one 
considers a horizontally stratified fluid under the action 
of gravity; then the question s What is the role of curvature 
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on the stability of the flow? does not have a satisfactory 
answer till today# In fact< the validity of the Rayleigh's 
result in the presence of stratification has not been 
established so far# 

As early as 19 33# Synge gave a generalization of 
Rayleigh's theorem. Later on# Yih (1974) gave his sufficient 
condition for stability^ which is the strongest known genera- 
lization of Rayleigh's theorem. Unfortunately# both these 
results are not as siirple as the Rayleigh's theorem. Trorpe 
(1969) has found an excinple of a basic flow which is stable 
in the homogeneous case by Pjortoft's theorem# but which is 
unstable in the stratified case# In this chapter# we give a 
siitple generalization of Hoiiand' s restilt to the stratifi^ 
case. 

Howard (1961) was the first to give an estimate for 
the growth rate of an unstable mode of a stratifi^ fluid# 

% 

For homogeneous fluid# his estimate reduces to Itoiland's 
estimate of the growth rate# Recently# Makov and Stepanyants 
(1984) gave a generalization of Howard's estimate taking in+p 
account the total depth of the fluid layer with velocity shear# 
All these estimates do not take into account the curvature 
effects on the gixjwth rate# 

For a homogeneous fluid# Sattinger (1967) gave an 
estimate of the growth rate of an unstable disturbance involving 
the curvatiire of the basic flow# This has been inproved by 



25 


Craik (1972) In this chapter* Sattingei/s estimate is 
generalized to the stratified case as it is not possible to 
do the same for Craik's estimate# 

Recently Kochar and Jain (19 79) proved that the complex 
wave velocity of unstable modes of a statically stratified 
fluid lies inside a certain semi-ellipse in the upper half- 
plane* whose major axis coincides with the diameter of Howard'^s 
semicircle* while its minor axis depends on the stratification 
parameter viz. the Richardson number* We show that the 
unstable modes lie also inside a semiellipse region whose 
major axis coincides with the diameter of Howard's semicircle* 
while its minor axis depends not only on the stratification 
but also on the curvat\are of the basic velocity profile* 

Lastly, we extend the results of the previous chapter 
to the stratified case. We derive a necessary condition for 
instability which simultaneously gives Miles' criterion* a 
range of c^^. and c^^ and also takes into account the oirvature 
effects on the stability of flows# It is interesting to 
note that these resxilts are inprovements over the results of 
Banerjee et al. (1974* 1978). , 

2. Holland' s Griterion . 

The Taylor— Goldstein system is 

+ ( 1 ) 

(U-c) ^ 



with the associated boiindary conditions 


H(y) = 0 at y = » (2^ 

If (l) is multiplied by (coitplex conjugate of H) and 
integrated over 9^/ after an integration by parts 

/ p( iH'i^ + + / IfHlil jh}2 - / _ 0 

(U-c)2 


(3) 

The limits of integration and the infinitesimal length dy are 
dropped throughout for convenience. The real part of (3) is 


/ p (iH^r+k"^iHr) + / 


julcp 


2 pN^[(u-cj^-c^3mi2 


mr-j 


lU-cl 


(4) 


Using the fact (U-c^)^ - c^ s= lU-cl^ - 2c^ (4 ) can be 

rewritten as 

(PU')'(U-c )-pN^ 2 7 7 9 •> 

/ : — IH}= - / p(iH^r + k^fHr)-2c^ s — — 

lu-cl ^ ru-Gi^ 

(5) 


The right hand side of (5) is negative* Therefore, a necessary 
condition for instability is that 


(pu^)' (U-c^) < pN^ (6) 

atleast at one point in iYifY 2 ^* Since for a homogeneous fluid 
2 

p s= constant, N =0, the condition becomes 
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atleast at one point in the flow domain* This is Holland's 

criterion for the Rayleigh problem* Thorpe (1969) has shown 

that a stable flow of a homogeneous fluid with velocity 

profile U(y) = c - sinh y becomes lonstable, if the density 

7 o 

profile is given by N ' = tanh' y^ > 0* It is easy 

to see that these profiles do not violate condition (6)* 


3 » An estimate for the growth rate * 

From (4)^ we have 

9 0 0 (PU')' (U-c_) „ 

/p(iH'r + k iHr) = - / oj - — ^ 

lU-cT' 

.2 r _ ^2 


[(U-c ) 2 .- 4 ] 

+ / — — I — — i- mr 
iu-cr 


(s) 


Taking into account that 


- CPU')' (U-Cj.) < l(pU'3'J lU-c^l 


< 1(PU')'} (U 

- max max 


U . ) 
min 


and (U*c^)^ - c| < lu-cl ^ 
we obta-in from (8) 


; PdH'iWiHi^) < HPa'yi I 


iHl 


max max min ju-cl^ 




mr 


max jy^^j2 


( 9 ) 


Using the well known Poincare inequality, we have 
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2 

/ PdH^I^ + > P . (-—2. + y^) / IHI^ (lO) 

Using (lO) and the fact w < in (9)/ we get the 

iu~cr 

following estimate : 


,, o I(PU')'l (U -U . ) + (pN^) 
■,.2_2 ^ max max man max 

k ^ 

( 1 ^ ^ 
min ^ , 2, „ i2 


(ll) 


For a homogeneous fluid (p = constant)^ this reduces to 
Sattinger's estimate 


■1,2 2 ^ *max ^^maac " ^min^ 
jc c^ < 


(1 + 


71- 




) 


( 12 ) 


From ( 12 ) we see that constant vorticity flows of a homogeneous 
fluid are stable. But we can not make such a conclusion from 
(ll)» Therefore constant vorticity flows, which are stable 
when the fluid is homog®ieous, may become unstable whoa the 


fluid is stratifi^. 

4 » A Saniellipse Theorem . 

If we intraiuce the transformation H = {U--c)f, then 
the Taylor-Goldstein system becomes 

[p(u-c)^p']' + p[n^ - k^(u-c)^] p s= o (13) 

and FCyj^) = O = FCy^) 


( 14 ) 
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If the equation (13) is multiplied by F* (cortplex conjugate 
of f) and integrated over then the real and imaginary 

parts together with certain ingenius manipulations introduced 
by Howard (1961) irrply the Inequality 

[ (c^. - + c^ - (^^)^*] / PQ + / pN^ }FI^ < O (15) 

where Q = IF^l^ + k^'lFl^, a < U(y) < b- (16) 

Howard drops the last term in (15) which is positive to 
establish the s^nicircle theorem* Kochar and Jain (1979) 
found a positive/ lower and relevant estlnate of the dropped 
term to get their semiellipse theorem* However, neither 
Howard's thasrem, nor its generalization by Kochar and Jain 
take into account the dependence of the size of the instability 
region on the curvature of the basic flow* Now, we find a 
new estimate for the last term in (IS) to get a new sexii— 
ellipse theoren Incorporatjhg the ctirvature effect 

H = (U‘-c)F. 

Differentiating this/ we get 

IH'} > i lU-cl iP*! lU'J iFf 1 
Therefore/ / PlH-'l^ >/ piu-cl^“ IP'I^ +/ plU'l^ |FI^ 

-2/PlU-c! lUM IFI IF' I (17) 

2 

Let = / PU'^ IFI^ ; a / p ju~cl^ Q, J^ s ' 
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where J(y) denotes the local Richardson nuniber* 

By Cauchy-Schwcirz inequality 

J PlU-cl IP! IP'i < [/ PU^^JFI^ / PlU-cl^'lF' 

< be (13) 

Addition of / pk^lHl^ to both sides of (1?) and use of (18) 
gives 

/ PCIH’'!^' + k^lHl^) > (B-E)^ (19) 

Using { 9 )/ we get 


(B-E)^ < [ 1(PU')'^^^ (b-a) + (PN^)^^^] / IFI^ 


(20) 


Now, since / IFl^ < I 
we obtain from (20), 


pU^2 jpj2 


P . . 

min min 


j— 4 provided ^ 0, 




p U* ^ 

min min 


5 


< 1 + A 


£21) 


l(pa^)'l (b-a) + CPN^L_ 
where ma^ y* ^ q, 

■ J, ' luxn ^ 


P 

min min 


Therefore, / pN^ iF|^ > j B^ 

^ ■ r\ 


J E" 

■ - 0 . ^ * 

{l-hA)“ 
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J c\ 

> -2 — K I PQ. 

“ (1+A)^ 

Thus we have established the following leimia4 
Lemma 1 - 

For an unstable mode# we have 


/ PN^IP]^ > 


J c;‘ ; PQ 
O X 


{1 + [ 


'‘‘'“'’'Lax 


( 22 ) 


max 


P . V*. 

min min 


provided O. 

Let us use the estimate given by the l«nma of the 
last term in the inequality (15) to get 


{ (Cj.-^)%cJ [l + 


■l(PU')'l„,^(b-a)+(PN2)„^^ . , 

£ max max jl/2^ 2 

P ^ ^ 

min min 


-3 


- (^)^} / PQ < O. 

This implies the following theorem. 

Theorem 1 . 

For flows U with # O/ the oomplex wave velocity c 

for an tanstable mode must lie inside the semiellipse 




l(pu')'l fb-a) + ,,, , 

£ . "MBX - ^ max jl/2j 2 


-] <= 


^mih ^mln 


( 23 ) 
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In particular, for a homogeneous fluid for which = 0, this 
semiellipse theorem reduces to Howard's semicircle theorem# 
Jain and Kochar (1983) and Makov and Stepanyants (1984) have 
given generalizations of the Kochar-Jain semiellipsa theorem# 
Following them, similar generalizations can be given for our 
semiellipse theorem also# However, they do not s^m to throw 
any additional light on the role of curvatiire on the stability 
of the flow and hence, they are not given here# 


5 » Reduction ■ and Unification Problem* 

Banerjee et al« (1974, 1978) have established a necessasry 
condition for instability which simultaneously gives Miles' 
criterion and a range of Cj, and for an unstable mode for a 
class of velocity and density distributions satisfying (i) a > O 
and (ii) (PU'')' 0 or > O everywhere in the flow 

domain# Gupta et al# (1982) have removed the restriction 
(i) and obtained the resiilts# All these results have the 
inherent drawback that they can not be aj^lied to honogeneous 
shear flows for in this case condition (H) in view of Rayleigh's 
inflexion point theorem implies tha t the system is stable# 

Now, we shall relax the very restrictive condition (ii) and 
thereby generalize the results of the previous chapter to the 
stratified case* 


The stability equation in terms of G = (U— c)^/^ F is 


[p(u-c)G']' ^ Pk^cu-c) t 


]G =0 
(24) 
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and the associated boundary conditions are 

G(y^) = O = 0(72^ C25) 

Theorem 2 * 

If (i) a > O and 

2 p 

(ii) f(y) = [(PU')' + — > O, 

(y2-yi)‘ 

for every y e 

then < X (c^ - ^26) 

2 P(iTH-l) (™ - , 

where X = max { tt — A — — } and m = - . 

Proof * 

Multiplying (24) by G* (conplex conjugate of G) and 
integrating over the range of y, we have upon integration by 
parts and making use of (25), 


/ P(U-c) [IG'I^ + k^'lGl^ ]+ I / (pu')' |G|^ 

P(2^-n2) , • 

“ ° 

Equating the real and imaginary parts of (27) to zero and 
cancelling Gj (> O), we get 


/ P (U-Cj,) [IG^I^ + k^lGl^] + I X (PU^)' iGl 


+ / 


p(U-c^)(S^-H2)IGr 

(U~Cj.)^‘ -I- 


= o. 


*(28) 
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I ^ t 1,2 1 ( 2 T . f ; ;_4 I /-I I 2 


; p [iG' r + k^jcr] + s 


p(n 2 -I5i) 


(U-G )^+C? 


iG\^ = 0 


(29) 


r X 


Multiplying (29) by rtK:^ and adding the resulting equation to 
(23)# we get 


/ P[U + (mr-l)c ] (IGM^ + k^'lGl^') + ^ S iPUn* IGI^' 


2 

+ S Q 

(U-c^)2 + cj 


(30) 


Since a > 0# and a < c„ < b (Miles# 1961)# it follows that 


[u + = b > 0 and [u-(iiH-l)cJ^ = < O 


(31) 


Equation (30) \jpon using (3l) and the Rayleigh-Ritz inequality 
gives 


2 

f(y)G? + 2P {u-(im-i)ci (2~ - , .2 

; _-i ^ ^ 4 IGl- < 0 (32) 

(U-C^)^ t 

Under the conditibns of the theorem# ( 32 ) clearly iitplies that 


^i - ^^“^r “ 

2P (lIH-l) 

where X = max { — 

[yi,y2] ^ 


N^) 


and this establishes the theorem. 


Theorem 3 * 


Uaier the conditions of th« 3 rem 2# if 


?^ < a+b) - [(^ atb)2 ^ (b~a)2]^2 


( 33 ) 


35 


then the parabola 


c; = X(c - — 

1 r rtH-1^ 


intersects Howard^ s semicircl« 


(34) 


Proof# 



2 


* =1 = 


(35) 


It Is easily seen that the parahola (34) touches Howard's 
semicircle (35)^ if 

^ = (^ ^b) 4 [(^ atb)2 _ (b-a)2]V2 ( 35 ) 

The value of given by (36) with the positive sign is 
rejected as it leads to c^ < a which violates a < c^ < b. 
Therefore 

Xe = a+b) - [(^ - (b-a)2]V2. 

Hence, If X < (^ a+b) - [(^ atb)=’ - (b-a)2]V2, 

then the parabola ( 34 ) certainly Intersects Pioward's semicircle 

(35). 

This establishes the theorem. 

Theory 4 . 

If (i) a. > O ar^d 

2 

(il) g(y) =[(PU')' - -. -- s ax ^ 
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for every y in [ 731 / 72 ]# 




2 P(in-l) C ^ 
where X = max { ^ 

[^1/72] lg( 7 ^i 


/2 


- N^) 


(37) 


and m = I > 1. 

P roof * 

Multiplying ( 29 ) by — mc^ adding the resulting equation 
to (28), and proceeding as in Theorem 2, we get the res\ilt* 


Theorem 5 * 

Under the conditions of Theorem 4, if 


bfa)^ - (b-a)^^^, (38) 

then the parabola 



t 



+ 


m-1^' 


(39) 


intersects Howard's semicircle (35). 

Proof * 

Follows by proceeding as in Theoaran 3* 

From Theorems (2) or ( 4 ), It follows that a nasessary 

j|2 4 

condition for instability is that J = — < 4 somewhere inside 

the flow domain (Miles, 1961), Further Theorems 2 - 5 in 
conjunction with Howard's semicircle theorem clearly show that 


if 
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a > 0 and either f(y) > 0 or g(y) < C, (40 ) 

for evary y in [7^/721/ then under conditions (33) or (38)^ 
Howard^ s semicircle bound for hetrogeneous shear flows can 
further be reduced. The inportance of condition (40) lies in 
the fact that under this condition the basic velocity pTOfile 
can change the sign of its curvature in the flow domain* 
Further^ for homogeneous shear flows they reduce to the results 
of Chapter 2. 

& • Concluding - Remarks * 

Cuirvature of the basic velocity profile plays an 
inportant role in the stability analysis of inviscid, inconpire— 
ssible plane parallel flows in the absence of stratification. 

It is through curvature -that a condition which is necessary 
as well as sufficient for the existence of xinstable modes is 
available. In the presence of stratification one won<ters 
whether curvature of the basic velocity profile has the same 
role to play? The present study has been motivated by the 
above q[uestion and all the results of this chapter explicitly 
incorporate the curvature effect on the stability of stratified 
shear flows* However, one still does not ]<uow whether the 
sinple Rayleigh^ s theorm is true or not in the presence of 
stratification. 



chapter - 4 


ON THE BOUNDEDNESS OF THE WAVE VELOCITY OF NON-SINGULAR 
neutral MODES OF THE TAYLOR-GOLDSTEIN PROBLEM 

3 • Int.roduct ion 

For a given flow and wave number, the normal modes 
may be divided into five classes, some of which may be eirpty 
(Drazin and Reid, 1931). The investigations of Howard (1961) 
and Kochar and Jain (1979) show that the conplex wave velocity 
c = Cy -f ic^ of unstable (c^^ > O) and their conjugate danped 
(Cj^ < O) modes of the Taylor“Goldstein problem lie in a boxmded 
region in the c^c^ -plane.' The wave velocity of singular 
neutral modes by their very definition (c^ = O, U-c = O for 
some y in [yj/ Yj] ^ lies between the minimum and maximum values 
a and b respectively of the streaming velocity U« However 
the definition of non— singular neutral modes (e^ - 0 , U-c O 
for any y in irrplies that their wave velocity is 

either bounded below by b or boxmded above by a* Therefore 
the possibility of their unboundedness still exists. In this 
chapter we rule out this possibility by proving that the wave 
velocity of an arbitrary non-singular neutral mode of the Taylor- 
Goldstein problem is boxinded* However, if the fluid is 
coirpressible, the wave velocity of an arbitrary non-singular 
neutral mode may be unboimded as shown in the next ch^ter. 
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2 # Establ ishment: of the Result » 

For a non-singular irode, ws have c = and U-c^ ^ O 
in illustration by Rayleigh's (I883) exanple, 

we take the Taylor-Goldstain equation under the Boussinesq 
approximation* Therefore^ the Taylor-Goldstein equation may 
be written as 


H» + =0 

Cu-o_)’- W - V •' 


(1) 


and the associated botindary conditions are 


H(yj) = 0 = H(y2) 


( 2 ) 


For a real eigenvalue c = the eigenfxinction may be takoi 
as a real valued fvinction* Therefore, multiplying equation (l) 
by the eigenfunction H and integrating the resulting equation 
with the help of (2)/ we have 

/ + IHl^ (3) 

Using the Rayleigh— Ritz inojuality (Shultz, 1973 ) in equation 
(3) and dropping the positive tem involving the wave number 
from the resulting inequality, we have 


/ 


[(U-c„)^ + U"(U-c_) - ]— Jiil— s 


< Q 



(4) 



where 
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Inequality (4) iiiplies that there exists e [y^, y^ ] such 
that 


f(c ) = - [ 2 U(y )a^ + V" (y^)] c, 


s r 


+ [aV(yg)+u(yg)u"(yg)“N 2 (yg)] <0 ( 6 ) 


2 

When N > 0^ the discriminant of f(c^) = 0 ^ namely/ 

\ 2 2 2 

0 (yg) = U*' (yg)+ 4 a N '(y^) is positive* Therefore/ it follows 


from inequality 

( 6 ) that 

\ 



A < 

Cj. < B 


(7) 

where A = UCy^) 

u"(y-) 

- Vb(y ) 


+ ® 

— -5 

( 8 ) 


2 a 


and B = U(y ) 

U*'(Ys). 

+ '^Cy_)’ 


+ 2— 

’ S 

2 

( 9 ) 


From these/ we have the following theorem* 


Theorem- 1 * The phase velocities c =5 c^ of non-singtilar 
neutral modes of the Taylor-<?oldstein problem are bound®! 
and 


* ^''1 rr ,-? I ! +^ 2^2 

- 2 — ’^C^N EH- „ 2 


2a 


2a 


( 10 ) 


2 ' 'ifr' ' 

rhere a =* ^ 


(Yj-yi) 
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It is interest. \ng to note that the bounds givenby (lO) do not 
depend, upon the wave number However^ if we do not drop 
the term multiplied by k" to arrive at the equation (4)/ then 
we get the following theorem* 

Theorem 2# 


For the phase velocities c = of non-singular neutral 


modes 


lu 


tt\ 


a + 


‘min 




2(a^ + k^) 


1 Cj. < b + 


lU' 


^max 


+ 4(a\k^)N3 


max 


maX' 


2Ca^ + k^) 


( 11 ) 


From equation (6 we see that a necessary condition for the 
existence of non-singular neutral inodes is that U"^+4a^^ > O 
atleast once in the flow domain* But, Miles (1961) has ^own 
that if N (y) < O throughout the flow domain then non— singular 
neutral modes do not exist* However, for homogeneous flows, 
our result inply that non-singular neutral modes do not exist 
if u(y) is linear* 

3# A Simpl e ^ II lu st ration » 

For U 2 O and P = Pq \d 1 e 3 re 0 > O is a constant, 

Rayleigh (I 883 ) showed that the non-singxalar neutral modes 
are given by 

2 

c = + [ Y -2 3^ / ° = 1^2, .. . 

^ k^ + n a 


( 12 ) 
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where 2sr^ - rro 

- positive constant, it fo]-, 

(12) that follows from equation 


•^ < c < N 

“ r < a 


(13) 


3nQ this is orecj cioi tr i. 

present case, equation (n) 


2 1 ; 1 
~ r N- - ^2 .,2 1 

r -^]2 < c ^<[_ n __ j 2 


a 2 +k 2 - 


a^+k^ 


(14) 

and we see that these bounds coincide with th 
Of C, given by («) „,en „ = 1 . 

^SS ClMinq Remar-Vo . 

We hQ.ve shown th;^+- 

neutraa rnirt « ™''s-velocity of non-singular 

neutral modes of the layiorK^oldsteln probl«, . . 

Hence, for the laylor-Kloldsteln p„bi«„ 

Of any nom« n..o Is hounded. furthlrT 

o^alned In this chapter are the hest::sti IT 

Illustrated, This bounda^ 

~^nxs Doiindedness of thf:^ 

singular neutral modes Is ^ ''moetty of non- 

only, as in the T ^-"Presslhle fluids 

in tne next chapter we show that they mav h . 

for conpresslble Shear flows. . ^ 



CHAPTER - 5 


stability of gasdynamic shear flows 
i* Introduction > 

Some general stability characteristics of gasdynamic 
shear flows have been obtained by Lees and Lin (1946)* Then^ 
Echart (196 3)/ Bl^Imen (1970) and Chimonas (1970) initiated 
studies on the general aspects of the stability of parallel 
shear flows of an inviscid fluid with particular reference to 
the effects of coitpressibility on the mechanism of instability# 
The extension of the Rayleigh stability criterion and ii>ward's 
semicircle theorem to conpressible flows/ obtained by Lees and 
Lin (1946) and Eckart (196 3) respectively/ have each bemi 
rederived by Blumen (1970). The Miles' th^ran and Howard's 
estimate for the growth rate of an mstable mode for stratified 
shear flows have been extended by Chimonas (1970) to stratified 
coitpressible shear flows* Dandapat and Gupta (1977) have 
studied the stability of magnctogasdynamic shear flows# 

In this chapter/ we follow the work of Blvimen (1970) 
and obtain general stability results# We have shown that shear 
free conpressible flows support supersonic waves but not subsonic 
waves# These supersonic wa-^s are non-singxilar neutral modes 
whose phase velocities do not lie inside any boimded region in 
the G^c^-plane# We have found an Instability region for subsonic 
disturbances vhich depend on the iiach number/ wave number and 
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depth of the fluid layer- For unbounded flows, we have 
rederived Hoiland^s estimate for the growth rate of an unstable 
mode and have also found a sufficient condition for stability 
to supersonic disturbances. 


In this chapter, we use the stability equations in 
their non-dimensional form so that the Mach number characterises 
the conpressibility. Also, we use the stability equation in 
terms of the perturbation pressure. Since, these were not 
introduced in the first chapter, we shall do it here. 


We shall consider the linear stability of the basic 
plane parallel flow U(y), of an ideal gas moving in the x'- 
direction, with transverse variations along the y'-axis. For 
sinplicity, the basic thermodynamic state is assvimed to be 
constant and is chciractsrized by the sound speed 


where p and p are the pressure and density re^ectively and ^ 
denotes the ratio of specific heats, Superpos«i on this basic 
state are small disturbances in the Cx'', 3 r^) coiiponsnts of 
velocity, (u^v'), and pressure p'* 

Non-dimensionalization is carried out by Introducing 
a velocity scale V and a length »3ale L, which are characteri- 
stic of the transverse variaticms of the basic current U, We 
then define the dimensionless coordinates, time, velocities 
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and pressure as 


(x^y) s (x^,y')/L , t s , 

j-i 

u S u/v, (u,v) S (u^v^)/V, 

71 = p'/pV^* 


“1 

! 

I 

I 

i 

I 

1 

1 


( 2 ) 


The Mach number is M H •— . Then the basic system of 
linearized inviscid equations becomes (e*g» Betchov and 


Criminal e, 1967) 



^t + Uu^ + V Uy 

= “ ^x 

(3) 

^t + % 

It 

t 

(4) 

+ u 71^) + 

"x + ''y = ° 

(5) 


where (3) and (4) are the momenttm equations and the ajuations 
of mass continiilty and entropy conservation have been combined 
to yield equation (5)» 


form 


Each wave disturbance will be represented in the 


q = q(y) exp [ik(x-ct)J (6) 

where q is u, v or jr, k ig the x-wave niimber and c = c^-f-ic^ is 
the conplex phase velocity# The stability problen is to 
determine the coitplex eigenvalues c lawier the conditions that 
the wave disturbance (6) satisfiesS the linear equations and 
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the boundary conaltions separately^ Instability corresponds 
to ^ ^ consequently an exponential growth of the wave 
disturbance at the rate of k c^# 

We shall make use of the differential equations for 
the aiTplitudes v and 7 i respectively of transverse velocity 
coiTponent and pressure « These equations inay be obtained 
from (3)/ (4)/ (5) and ( 6 ) in the form (after dropping the 
C arret ) 


(U^c)y’' - U'^y 

and 

(U-c)7l" “ 2^'7i' + ]C^‘(U-C) [ 1 - M^'(U-c)^] 71 = 0 (S) 

where primes denote differentiation w.r»t» y* 

Following Lees and Lin (1946 ), the disturbances 
associated with (7) or (8) can be classified as subsonic, sonic 
or supersonic depending on whether the relative phase ^eed 
lU~cl is less than, equal to or greater than the inverse of 
the Mach number 

i,e, lij~c! I (9) 

The subsonic disturbances are the counterpart of the so’-called 
inertial modes, which are soltitions of (7) in the limiting 
case M = 0 (the Rayleigh stability equation)# The physical 
significance of the sonic distua*ances is apparently not clear 


- k (U*-c)v = 0 


(7) 
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and «111 not be coneldered further. Plnally, auperscnio 
disturbances correspond to compression or sound waves. 

On rigid boundaries the normal velocity must vanish 

Thus, from (4), we have 


V = 71' = 0 (y = y^,y^) 


If the fluid is unbounded, then the boundary conditions 

7I=V=0 (y-* + oo) 

^ * .^i^^sonic and Supersonic Waves* 

In the absence of basic shear (U 5 0 ) (7) Income 


(10) 
become 

( 11 ) 


1-M C 


0 


( 12 ) 


Now, from Howard's semicircle theoran (Bltuuen, 1970 ), it 
seen that U = 0 is stable* Therefore (12) with (lO) can have 
only real eigenvalues* 


Lemma- 


In the absence of basic ^ear, subsonic waves do 


not exist* 


Proofs 


Mxiitiplying (12) by v and integrating over {yj,y ) 
gives after an integration by parts using (lO) 


,2 

( “C ) / .1 . 


2 2 

HI* ‘jNip IS* 

l-M^? 


(13) 



The limits of integration 
are dropped throughout for 


and the infinitesimal length dy 
convenience# 


It we consider only subsonic waves, then > o 

and we see from (13) that c = 0. 


This proves the lemma. 

Lemma 2 # 

absence of basic shear, supersonic waves do 

exist. 


Proof. 


as 


For c ^ 0 and 1 -M^c^ < 0, 


equation (12) may be written 


+ k^(M^c^^l)v = 0 
Now the function 

V = sin [ k (Vc^-i) (ywy^) ] 


( 14 ) 


( 15 ) 


satisfy the equation (14) and also the boundary condition (lO) 
whenever 


r - T- g ■ 

M ^ C ‘^.1 


- n7i ' 


(n = 1 


O 

J? JF 9 0 


• ) 


That is, whenever 
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This proves the lemma* 

Combining lemmas 1 — 2 , we have the following theorem* 

Theorem- 1 * . 

Shear free basic flow (U — 0 ) supports supersonic 
waves, but do not support subsonic v/aves. 

Remark 1 * 

Rayleigh (I 883 ) has shown that shear free basic 
flow supports internal gravity waves (in the presence of 
gravity). Here, we have shown that compressibility also 

new waves in shear free flows and that these waves 
are supersonic waves* 

Remark 

For inconpressible shear flows, we showed in the 
last chapter that the wave velocities of non-singular neiitral 
modes are bounded* But, in the presence of compressibility, 
this result is not true* An exanple is given by the supersonic 
waves discussed in lemma 2, which are non-singular neutral 
modes and their wave velocities given by (16) are unbounded* 

4* Instability Region for Subsonic Disturbances * 

Eckart ( 1963 ) and Blumen (1970) have shown that the 
coxiplex wave velocity c for any unstable mode lie inside 
certain semicircle in the tpper half of the c^Cj|^-plane* But 
their semicircle does not involve the Mach number and hence 
is the same as that of incompressible shear flows* Considering 
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only subsonic disturbs-nceSf we iirprove upon this result by 
finda.ng an instability region which depends on the Mach number 
and lies inside the semicircle of Eclcazrt (1963)# Interestingly^ 
this new instability region depends also on the wave number 
and depth of the fluid layer# 

T heorem 2 # 

The corrplex wave velocity c for any unstable subsonic 
mode must lie inside the semicircle type region 




M 


^31^(1 - I M^(b-a)^)^.c^ 


4 [3lW(y2-yi)^(l -| M^Cb-a)^)^*] 


< (^)2 


(17) 


Proofs 


For c. > 0^ the transforrration 


V ^ (U-c)P (18) 

reduces the original system (7) and (lo) to the syst®a 


.(U-c)P' O _ 7 

- k^(u-c)'‘F = 0 

1 M^(^c)^ 


(19) 


(20) 


and F = 0 at y = Yi^*Y 2 

Following the standard procedure (Howard^ 1961), we can obtain 
the inequality 
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[(Cr + cl - Q 


ll-M^(U-c)^I^ 


Where Q L ^ . 1 . 1 , , + k^lPl^, 


< 0 


a < U(y) < b. 


( 21 ) 

( 22 ) 


Dropping the last term which is positive, one can obtain the 
semicircle theorem# But, that is precisely the term which 
involves the coirpressibility effects explicitly# So, we shall 
find a lower and relevant bound for the ternv though only for 
subsonic disturbances# 

For subsonic disturbances, we have MlU-cI <1. 

Therefore ll^^(U-c)^I < 1+M^lU-cl^ < 2* 


Therefore j , ,^,,2 

ii-^^(u-c)^r ” ^ 


(23) 


From the semicircle theorem, we have, for Cj^ > O, 


lU-cl 


(U-Cj,)^ + cj < (b-a)^ + a I (b-a) 


Therefore, for subsonic inodes 


!l-M^(U-c)^l > Il^^'l^cl^f 


I ! 1"^ " 


,106208 


>1 - f ■ 
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Therefore, / Q < / £. 


ip'r 


[l - I 

Since by the seinicircle theorem 




(24) 


(Or - + oj . 


2 , 

< Ojr we have after using (23) 


and (24) 


[(c„ - 


J!j^) 2 , 2 f b«»a A 2 I . p 

-y-; + ] / [ 


IP' I 


(1 - 


2 A + k2|P!2] 


2 

+ — ^ / iP'r < 0 


(^s) 


Now, by using the Rayleigh-Ritz inequality 
/ IPP < — S-ji — / 


n 


in (25), we get after dropping / IP'I^ 






(1-| M^Cb-a)^)^ 


2 n- ] 




M^cf 

+ < o* 


Equation (1?) follows from this* 


(26) 


For inconpressible fluids M s o and the above result 
reduces to Howard's semicircle theorem* Furthermore, if either 


-* oo or (y 2 *’yi^ * “> then ^sof the alxjve region reduces to 



the semicircle. Bat, for finite values of and (y -y )2 

3 ilc 

and for non-zero, M, the instability region gets raluced# 

5 ■’ Stabilitv Analysis For Unboxmded Plnws . 

If Cj^ 0, and U is finite, then (S) may be divided 

4mm O 

by (U-c) ,with the result that 

i 

I 

[(U— c) 71 ^]^ — - k^ [(U-c) ^ ] TT = 0 (27) I 

Multiplication of (27) by ji* and integration over (- 00 , 00 ) and 

■ I 

application of (ll) yields j 

; [(U-c)"*^ Iti' + k^ { (tJ-c)“'^-M^} 0 = 0 (28) 

Following the standard procedvire (Howard, 1961) we obtain 

[(c^ - ■^)^' + cj - (^)^3/ Q + k^M^ / < 0 (29) 

where Q =* lU-cl ^ ( Iti'I ^ + k^Iirl^) ( 30 ) 

Since k M* > 0 and Q > 0, (29) gives the soraicircle theorsn* 

As noted by Blumen ( 1970 ), we see from (29) that 
increasing values of (kM) plays the sar® role in reducing the 
allowable range of unstable la^des as increasing values of the 
Brunt— Vaisalla frequency in Howard^ s result# 

Since we consider only unstable inodes, the transfonnatloni 
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±s «ell defined. Use of this transfonration to ( 8 ) results 
in the following equation in G : 

[(U-c)-l G']' + I [(5-o)-2ij']' G - I u' 2 (5 .c)-3g 

+ k2(u-o)“2 (3 ^ k2M2(u-o)G = 0 (32) 

and ( 11 ) gives G(y) = 0 (y ^ + oo) , ^^ 3 ) 

If ( 32 ) is multiplied by G and integrated over the 
flow domain, we gat after Integration by parts idienever necessary. 


S (U-c)"^ [IG'I WlGl^] + I ; (U-c)“2 u^GG*)- 


/r: ^-3 


+ / (U-c)*"'^ iGi^ - / A^(U-c) IGI^ = 


( 34 ) 


The irnagiii3.ry part of this gives 


(U-c„)U^(GG*)^ U'^(U-C^}2 IGI^ 

^ + * — :z:—r— + 1 


I U-C i 


lu-cl 


lu-cr 


, r 2 „ 2,^,2 , ^ 

+ S + / k M IGI - / ; y S! 0 

lU— cl ‘ 4 lU— Cl 


( 35 ) 


Considering the fir^t integral, we see that 

. j j 2 ( y-c^ )U^ ( GG*) ' y' ^ . 1 ^ 

lU-cl lU-cl |U~cl 


* 1--Z — : :■= — — — j 


lU-cl 


lu-cl 
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^2IU'| lU-c^l IGI IG'I (U-c^) U' (GG*)' 

L ^ -z: — = 7 — + ~ .. ' m ' 

lu-cl 


lU-cl 


]. 


(36) 


Concerning the second square bracjcet, we observe that 


(U-C ) U^(GG*)' lU-cJ lU^l I(GG^)/| 
r 5 < — 


lu-cl 


lu-cl 


2:1U-Cj,l fu^l IGI IgM 
lU-cl^ 


(37) 


Hence the expression on the left hand side of (36) is a non- 
negative quantity p (say), so that (35) can be rewritten as 


/ 



+ 


— sy) + / (yr yr - ■ ) lor 

lu-cl'^ 4iu-cr 


s 0 


(38) 


This iitplies that 


S 


k^lGl? ^ . 
lu-cr “* 



4iu-cr 


(39) 


From this,^ we get the estimate for the growth rate of any 
unstable disturbance as 


This estimate was first obtained by Holland (1951) for 
incoxtpiressible homogeneous shear flows and can be obtained 
as a special case, from (4*3) of Chimonas {1970)» 
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Equation (38) also gives the following theorem- 


Theorem 3< 


A 


necessary condition for instability is that 


2 2 T*T^ ^ 


4lU-cl'^ 


(41) 


atleast once in the flow domain < 


This result essentially means only that the flow t!(y) 
is stable to short waves. 

Now, for supersonic disturbances, we have !U-c)l > 
Hence, for supersonic disturbances, (38) inplies that 


2 2 

S < S (4 5'“ - k^"') , 


TTf ^ 


(42) 


From this* we get the following theorems* 

Theorem 4 * A sufficient condition for stability of the basic 
flow U(y) to supersonic disturbances is that 


2—2 ? 


(43) 


throughout the flow domain* 


Theorem- 5 * 

For an iinstable supersonic disturbance, an estimate 
for the growth rate is gi'fei. by ; 
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rr#2 


c" < (S^ - m-2) 


max 


(44) 


The estimate (44) is an irrprovement over the estimate ( 40 ), 

6» Goncludinq Remarks > 

We have shown that shear free basic flow (u s 0 ) 
support supersonic waves but do not support any subsonic wave* 
These supersonic waves are non-singular neutral modes whose 
wave velocities are not bounded# For subsonic disturbances the 
instability region given by the semicircle. theorem is further 
improved# The new instability region lies inside the semicircle 
and depends on the Mach number^ v/ave number and depth of the 
fluid layer# For unbounded flows, we have rederived Holland's 
estimate for growth rate of any unstable mode. Furthermore, a 
sufficient condition for stability to supersonic disttirbances 
and an inproved estimate for the growth rate of an unstable 
supersonic mode are edso given. 

The sufficient condition for stability to supersonic 
disturbances of an unbounded shear flow is similar to Miles' 
(l961) theorem for stratified inconpressible shear flows# Also, 
as ronarked earlier, the existence of supersonic waves in shear 
free basic flow is similar to the existence of internal gravity 
waves# These results seem to luply an analogy between 
coirpressibility and density stratififfatdon under gravity# 



CHAPTER ^ 6 


stability op stratified compressible shear flows 
1* Introduction . 

In thxs chapter, we study the linear stability of 
compressible shear flows in the presence of gravity. The 
stability equations have already been derived in the first 
chapter# Due to the mathematical coirplexity of this problem, 
no systematic study of this problem has been made so far. 

Eckart ( 1963 ) extended Howard's semicircle theorem to 
stratified corrpressible shear flows. Later, Chimonas {l9TQ) 
extended Miles' theorem and Howard's estimate for the growth 
rate of an unstable mode, to compressible flows* These results S 
were not simply found and needed considerable ingen it y. 

For incorpressible stratified shear flows, Kochar and. 

Jain (1979) have shown that the coipl ex 'mve velocities of ! 

■ ■ ■ 

unstable modes lie inside a certain setniellipse, rather than 

■ .. ■.'■■■ . ■ ' ■ 

Howard's semicircle, vdiose minor axis depends upon stratification 

■ 

This result has further be^ generalized by Jain and Kochar (1983 
and Makov and Stepanyants (1984)# Corresponding results do not 
exist for compressible flows* Furthermore, nothing is known 
regarding the role of the basic velocity profile 

on the stability of stjEafc|fied oorpressible shear flows. In this 
Chapter,, we 
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First, we TODdify Chimonas^s proof of theorem# 

This enables us to inprove upon the instability region given 
by the semicircle theorem for subsonic disturbances and for a 
class of supersonic disturbances# The instability region for 
subsonic disturbances, which is given in this chapter, depends 
not only on stratification but also on the wave number and the 
depth of the fluid layer# Furthermore, this region reduces to 
the line c^ = 0 when u;;^ / 0 and | - in 
theorem# 

To study the curvature effects on stability of stratlfie 
corrpressible shear flows, we make use of an approximation due 
to Shivamoggi (1977)* Under this approxiu^tion# we assx-ne that 

« 1 and « 1 so that their product can be neglected in 

^ , ^ 4 - •Ions# Under this 

comparison to unity, Iri the stability equatio 

g found an instabllit 

approximation and for subsonic modes we have 

1 involving the 

criterion and estimates for growth rate, an 

Also, we could extenc 

curvature of the basic velocity profiles# 

theory# 

many standard results of incor^ressible fio 



Modified Proof of Miles^ Theory # 

^ js g 

The stability equation in terms o£ 
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(2) 


and the boundary conditions are 

G(yj) = 0 = GCy^) 

For an unstable mode, the imaginary part of the equation, 

obtained by multiplying (l) by and integrating it using (2), 
gives 


II - 


Ig'|2 + 1 ig|2 




% 


U '2 , g ,2 


T 

f*. 


>. (GG*)' 

^ ^ * _ ' ' '*1 


. 2 
4!^ 


t / (n2 . - 


,2 r IGI 


= 0 


lU-cl 


(3) 


At this stage, we deviate from the analysis of Chimonas (1970) 
and rewrite equation ( 3 ) as 


(U-c^? 


11 


^ c(l f (IG'I^ t 


IGI^ 


■) 


4a: 


+ -4 (JG'l^ + ■— T- ~ ) - 




7 


4a; 


* 



(4) 


(5) 
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Therefore, 

(U-c U'^IGl^ U'(u-c ) (GG*)' 

(1 + (IG'I^ + — ;r-) - - — -I— — - 


4a' 




r 1/2 (u-c i ja,j 13, 

= [{1 + — iG'i - fi + — —2r— 3 


"2 ^ ^ • - I A T 2 ~ 

% % 


(U-n)^ lU'l IGI IG'I U'(U-c_) (GG*)' 

+ (1 + -— ■ £ - -l - ——4— — - 


’••Jt 




> 0, (by 5))* 

Therefore, for an lonstabie mode, (4) inplles that 

U'2|G|2 _ 


II . lMr§I-l2 


”7^) + ^ 
4a* 


' Y ' "" ' (IG'I'* + — T ^ «A j 


I }Gl^ 


( 6 ) 


^ IU-kjI 
2 

This is impossible if everywhere in [yj[»y23* 

This proves the following (Miles') theorem# 

Theorem 1 » • 

A sufficient conclitlon for stebllity is that N= > i n'^ 

throughout [yj#Y2^* 

3. Instability Regions # ; 

. ’ , '■* .V * ♦ 

. - ■■ ■ 




following theorem# 


Theorem- 2# 


Now, wen 

■' '-"'SS 

T’ fer" ^y unstatae supers>nic 
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mode satisfying ^ lie inside the senlelllpse region. 

in the upperhalf plane given by 


( ^ ^ 5t+b \ 2 , 2 r . 

(c^ “ [1 + 


4J, 


[!+(! + 


] 


% 


min 


< 


( 7 ) 


where a < U(y) < b* 


Proofs 


The stability equation in terms of P is (Comes from (23)^ 
Chapter 1 ) 


[ 


r(U-c)^ F' 


. (U-c)^ 

1 •• 


rk^(U-c)^ F + rN^P = 


0 


Multiplying this equation by P* and integrating it over (y^^yj) 
and following the standard procedure, we get the inequality 


[(Cr - / rQ + J 


b»“a\2 ■ 


IPl 


r lu-cl^ IP'j^ 


+ I 


2 e . ' (¥- e )?,2 “ 

a* 11 ^ ' T” * 


< 0 


( 8 ) 


where Q = ^^4Fi\ a < 0(y) < b. 




11 - 


'2, f-;i- 




we have' tixm (6 ) 

a# 



4 c|.; r Q. 

[l + (1 + 

^ %min 


( 12 ) 

I 


Using the estimate (12) in (8), we get (?)» 

This proves the theorem# 

Theorem- 3 * 

The corrplex wave velocity c for any unstable sxibsonic 
mode must lie inside the semiellipse type region# in the upper 
half plane given by 




(13) 


< ^ 

^ 2 ^min -min^-'- - 
where a « ^ 7 7 


Proof 4 


Sinceir^ for -if# 
,2 2 


f’-iHSjpfcable modes# 




Theirefore# 


for an 
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unstable subsonic modej, we have from (6) 


c2 


II 






^ i - J.) J 


r y'2 JGI^ 


° lu-cl^ 


o r c . I U-C! I 7 

Let Tc^l 1 , 1^1. 1 2 ,pj; 


2 , _ (U-c)^2 

1^,11 


, rc? lU'l^ IPI^ 

and B, « / :_-± 


11 - 


(U-c)?,2 


r lU-cl 


Then^i proceeding as in the proof of the previous theoareia. 


E 


2 

2 ^ ^ 




Use of Rayleigh -Ritz inequality shows that 


r C^. IG'!^ 


ajfr 11 - — 

—a 

7 


(U-c) ’ ,2 


a*. 


> y 


c: riGl 


^ 2 

?*min ^min 
where K = — — 


I M *: 5 


Use of (17) in 1 


^axf ^ • ' ' 


(14) 

► 

•} 

we get 

(15) 

(16) 

(17) 
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- Jo' . 


(18) 


r IP'I^ 


If ^ S 


II 


/ i r 1P| / then use of (iS) gives 


-i-n •? 9 ^ ) (^ '*• J ) 2u U' 

/ - S9 . '^ = + k^' < ; . / i. : ,. :? . : .. :. , ; .., : , P' -^max- - it 


2 

•max 


r IPI 


/,2 , , 2^ 4 

(A + k ; c. 


(19) 


We have/ » 


r |y-clQ 

— r~- ” - 


r.c^. . 

a 2 u-i 2 ::£l!| 2 |u~ 

a* 


(V^ + k^) 


u' ^ _ 

max max 


(1 


5 (b~a)M 2 

*. i i. M. ..— ; 


4 a;^ 


min 


which can be rewritten as 


U'?- 


max 


Cv\k^) ~ 


~TT 2 

^i ^min 


S(b-a)M 2 

7 T 2 ~ e 2 
^^^min - -i 

l 2 • 


b; 


Using the inequality (IS)/ we get 



( 20 ) 
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Using this in (19), we have 


_6 


/ rQ 


/ rlPl 


f 1 . S(b-a)'-\2 rl . /I ' _o 4 

^max <1 - “a > [j + 2 /I - C| - J„) 


°i 4min 


(21) 


Therefore, / rN^lPl^ > j u'? / rlPl^ 

- o min ' 

>. ^min ^ ^ 


4.ax <1 - Jo) rf + 2 /i - Jo f 


( 22 ) 


Use of this in (8), gives (13)# 
This proves the theorem# 


In these two theorems, we have reduced the instability 
region given by Howard's srafdoircle* The reduced regions dqpend 
on stjcetification through the miniiraim Richardson momber J^# 

The second res\ilt for subsonic distxirbances incorporate, not 


only the stratification but also the wave number and the depth 
of the fluid layer# Idrthermore^ when / 0# >^0 •* f " 
inplies -♦ o'*’ in accord with Miles' theory# 



which throws light on the 
feh^.Tocity profile on the stability 
However, if we maOce 
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Illustrating the role of curvature on the stability of stratlfie* 
conpressible shear flows* 


The stability equation is 


rH' 


1 - 


(U-c)^ 




rU<- - ■■ T/ 


1 ^ 




V lu-c)' 




(23) 


and the boundary conditions are 


H(yj) =0 = H(y2) 


(24) 


Following Shivamoggi (1977), we shall assume that 
^ ^i that their product can be neglected in 

3 l 

corrparison to unity* Under this i^proxinfiLtion, the stability 
equation (23) becomes 


[- 


rH^ 


ru 


'-r^Y Tn=^ - rA + = o 




<u-c)2 


1 - 


r- 


T 


1 - 


(25) 




Now, we shall prove the following theorems for subsonic 

a 


- <>yr9x>' . „ 
disturbances, for which 1 — > 0* 


Theoram 4 * 

A necessary condition for instability to subsonic 
disturbances is 








(26) 
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atleast once in ( 7 ^/ 72 )* 
Proof* 


Multiplying (25) by H* and Integrating it over (y^,y ) 


using (24), we get 
r IK' 


^ + 1 rt^lHl^ e j[— S!' mL 

. (u-c )2 ■' 

T 


U~—^] 




1 - 


/ ssdji^ ,0 

(U^)^ 


* 


(27) 


For an unstable subsonic mode, the imaginary part of it gives 


r,r .-r«^-..,, . -2«' !“fr.>.. IHI?. 


(V^J‘ 

1 


(u*-c^) + lu-cl^ 


(28) 


% 


The theorem follows from this# 

It is interesting to note that ss o, this resxilt 

reduces to Synge's criterion <1933)* Now, letting 


L = i[. 


* we have the following theorem* 


1 ~ 


Theorem- 5 * 

A necessary 
disturbances is 


In^abillty to subsonic 
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C j < ) 

i max 


(29) 


whenr^ver L has a definite sign over (y_y ) 

i 2 

Proof# 


Equation ( 28) can be rewritten as 
/ C[l(U~g^)-n]^' + (L^' c? - N^)} ^ 


0 


iU“cr 


(30) 


This irtplios that a necessary condition for instability is that 


2 2 2 
L C^ - < 0 


atleast once in the flow domain, 


(31) 


The theorem follows from this* Pof incoirpressible 
stratified shear flows/ the corre^onding result was obtained 
by Synge (19 33) and later, independently, by 'Ylh (1959). 

Theorem 6 » 

An estimate for the gro^'ddi rate of an unstable subsonic 
mode is given by 

rU^ 


J(- 


(b-a) + (l^^^rnav 

(U-^ r 


,2 2 ^ 

k o^< 




'mm 


(1 + 


‘S’ 


- 3 ) . .. 


Proo^ 


The real part of v;^. 
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[1 - ] 

d X* 


(U-c^) IHI^ 


iV-C J^ fU-cl^ 


+ / — — -^ ■ •' 

lU-cl^ 


7?^]' ^ 




lu-c| 


This gives 

^min^ (IHM WlHl^) < j {I( - ~ ~ - ^ . . , )M (b~a)-i-(riJ^ j^Mll 

(U-.cp2 max lU-cl' 

— i. 

a 

C33) 

Use of the Rayleigh-Ritz inequality in (33) gives the estinate 
(32). 


All the three theorems obtain^ using the approxiiration 
o f Shivamoggi involve the curvature of the basic flow profile. 

5* Instabll itv Conditions . 

In addition to the above resxilts involving the 
curvature of the basic velocity profile, we can obtain results 
generalizing many standard results of inconpressible flow 


theory. Now, we ^aUti;! 


.1^ semiellipse and 

^p^iiUl;^,ip«^ar Jain (1979, 1983) 


generalized sendellip^ :;t®P<^:i^’»; 9 |l_*C« 2 har axw 

^ ^ - •U • ■a74:7M/ 7: ^ ^ ■ 
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and the bound for the g,:owth rate of Matov and St^an^^ants (1984) 

The tranaformatlon H = (u^)p, „hen uaed in (25). leads to 
the equation 


r(U— c)"^ F' 


[- 


-^] *' rk^(u-c)^P + rN^P = 


(U-c^) 
1 ^ 

■ 


o 


and the boundary conditions in terms of P are 
P(y^) = 0 = PCy^) 


(34) 


(35) 


Following Howard (1961 ), it is easy to show that a necessary 
condition for instability is 


[(Cj. - •^)^* + c| * (^)^] ; rQ + ; rN^ < 


0 


(36) 


where 


Q 


I F^ [ ^ ^ 

(U-c_)\ 

[1 - — 


+ k^IPl^ 


(37) 


Dropping the last term, vrt:iich is positive, from equaticm (:^) 
gives us the semicircle theorem* Finding a lower and relevant 
estimate of the last teimi in (36) will lead to the semiellipse 
theorem* Before proceeding to do that, we shall prove the 
following theorem# 

Theorem- 7.* 



disturbances is that 



73 


4 [l - 


*mln 


atleast at one point in (y^^Yj) 


Proof. 


D^se of the transformation G = (u-c)^/^ 
to the equation 


I" in (34) leads 


r(U-c)G' 

r - 1 

•f , X ’ " 


(U-C_) 


d' G - 


rU^^G • ■ 
(u-c_)^ 


4 [1 - 5 —] (u^l 


- rk^(U-c)G + = 0 


and the associated jbotindary conditions are 


G(yj) = 0 = GCy^) (40) 

Multiplying (39) by G* (coitplex conjugate of g) and integrating 
it over ( 7 ^/ 72 ^ using (40 ), we get 


S r(u-c) t- 

[1 


j- + k^roi^) + 4 s , ] * iGi' 


I ^ 1 /-»f 2 


IG( = , 


a^: 


The imaginary part 


(40) 



/ [ 


[l ^ ] 

a 


+ rk^lGl ^ , m S , 


U 


,2 


(U-c^)™ 

4 [ 1 - —'T* ] 


, 2 l ri_i: 

'' j ^ 

i -5 3 


This is inpossible if N > 


U 


.#2 


(U-c^ 

4 [1 - — 


pr In'h.l 


Hence a sufficient condition for stability to subsonic 
disturbances is 


2 

N > 


U 


f 2 


4 [1 ~ 


(u-c^)' 


C 42 


throughout the flow domain * 

By the senicircle theorenv for an unstable subsonic wsKi 
a < c < bi» But, then 


(U - c) 


* ^%iin 


(43l 


Therefore, 


u 


^2 


U 


,2 


4 [1 - 


^■^in % 


Hence, a necessary condition for instability is 


C 44 )' 


< 


U 


* 2 


(U-C^) 


2 ~ 


U 


/2 


4 [1 - — ^ 4 [1 - 

^mln 
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atleast once In * 

This proves the theorem. 

The above theorem is true only for subsonic modes and 
that too -under an ^proximation » Furthermore, it is weaker 
than Miles'' theorem proved by Ghimonas (1970)» However, -the 
method of proof is the standard one (Howard, 1961 ) and does 
not invole any ingenius manipulations* So, it helps us to 
prove -the following theorems* 

Theorem' 8» 


The coirplex wave velocity c fpr any ^anstable subsonic 
mode must lie inside the semiellipse region, in the upperhalf 



< 


( 45 ) 


Theorem 9 . 


The coitplex wave velocity c for any unstable subsonic 
mode must lie inside the semiellipss-t-ype region, in the 
upperhalf plane, given by 


(c^ - 


a+bs 2 
2 ^ 


+ Cji^ + 


Jo Jjln 4 


u. 


/4 

'max 


[■ 'J ^ JoTi- 

4Cl.i|3h, ^ 

^*min 




Cb-a)- 

2 

^min 




^here — 


min 


^max 


(47) 



An estimate for the growth rate of an unstable subsonic 
mode is given by 


■ 2 

'max 


U'^t. [ 


Al< 


4 [1 

- ^^min 


- J 1 
a J 


[14. V - 

■'f 


I] 


9 reduce to those of Kochar and Jail 


1983) ^d 10 reduce to that of Makov and 

St^anyants (1984)* 

J'or the proofs of these three theorems# one has to 
ceed generally along the lines of the proofs of theorans 2 
3, and hence are not repeated here* 

^ ' ^-^P^^ludinq Remarks . ’ ^ 

In this chapter, we have studied the linear stability 
corrpressible shear flows in the presence of gravity. Due 
o -Che extreme mathematical difficulty, not much work has been 
'^^^^ne in this direction* Following Eckart (1963) and Chimonas 
<1970), we have attenpted to generalize many of the standard 
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results of stretlfled in compressible flow theory to the 
cortpressible flow case# First of all^ we have modified Chimoi 
proof of Miles' theorem,. This enables us to reduce the 
instability regions for subsonic disturbances and for a class 
of supersonic disturbances# The reduced instability region 
for subsonic disturbances d^end also on the wave number and 


depth of the fluid layer* It is interesting to note that 
this region reduces to the line c^ =: 0 as ** 4 *"/ when 


U 


r 

min 


/ 0 * 


This is in accord with Miles' theorem and thus Mile 


theorem has been linked to the instability region* 


The role of curvature of the basic velocity profile on 
the stability of stratified coir^ressible shear flows is also 
,?«]»r,-^yhii% 4 we^use an approximation due to Shivamoggi 
f.v: ''■'.{i97‘7)''*| Under .'^^roximat ion, we have generaiiz-ed ;an ' 

■ ' 4,' -.ft’’ ' ■ ■ 

instability criterion and an estimate for growth rate, 'v^ich 
were given by Synge (1933) for inconpressibl e shear flows. 

Also, the estimate given in theor^i 6 is a generalization of 
the estimate given by equation 12 and chapter - 3* In additions 

■ ' ’ .... 'I 


many standard results of stratified inconpressibl e flow theory 
have also been generalized. ’ 


However, many questions still ranain open. For example, 
the existence of waves in shear free stratified oorrpressifcle 
flows has not been proved. Also, the boundedness or unboiandedn< 
of the phase velocity of non-singular noatral jnodesl has not be< 
proved yet. 
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